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Kinematically forbidden processes may be allowed in the presence of external gravitational fields. 
These can be taken into account by introducing generalized particle momenta. The corresponding 
transition probabilities can then be calculated to all orders in the metric deviation from the field- 
free expressions by simply replacing the particle momenta with their generalized counterparts. The 
procedure applies to particles of any spin and to any gravitational fields. Transition probabilities, 
emission power and spectra are, to leading order, linear in the metric deviation. It is also shown 
how a small dissipation term in the particle wave equations can trigger a strong back-reaction 
that introduces resonances in the radiative process and deeply affects the resulting gravitational 
background. 



Processes in which massive, on-shell particles emit a photon according to Fig[T] are examples of kinematically 

^ ' forbidden transitions that remain so unless the dispersion relations of at least one of the particles involved are altered. 

' This possibility presents itself when particles travel in a medium or in an external gravitational field. 

C^^ '■ The action of gravitational fields on a particle's dispersion relations can be studied by solving the respective covariant 

^^ . wave equation. This can be done exactly to first order in the metric deviation 7^^ = (^^^ — 77^^, where rj^y is the 

r— 1' Minkowski metric .1-8]. The solutions contain the gravitational field in a phase operator that alters in effect a particle 

^ four-momentum by acting on the wave function of the field-free equations. This result applies equally well to fermions 

I and bosons and can be extended to all orders in ^^y . The calculation of even the most elementary Feynman diagrams 

^ . does therefore require an appropriate treatment when gravitational fields are present. While the inclusion of external 

i__j' electromagnetic fields has met with success in the case of static (Coulomb) fields [9% and can be easily carried out for 
static (Newtonian) fields, no systematic attempts have been made for relativistic gravity. The procedure developed 

'"^ ' below is intended to fill in part this gap and applies to weak, static and non-static gravitational fields. 

^ , Let us assume, for simplicity, that P in FigH] is an incoming fermion and that the photon i and outgoing fermion 

j«^ ' p' are produced on-shell. The solution of the covariant Dirac equation, exact to 0(7^1^), is [5| 

^. : ^{x) = -—{-l-|^'{x)V^-v^)e-'''-^^{x) = f^^, (1) 

1^ . 27n 

^^ . where I?^ = V^ -\-iV^{x), V^ is the covariant derivative, T ^{x) the spin connection and the matrices ^^{x) satisfy the 

^^ ' relations {7^(0;), 7" (x)} = 2gf^'^{x). Both r^(a;) and 7^(2;) can be obtained from the usual constant Dirac matrices 



by using the vierbein fields e^ and the relations 



rix) - e^^{x)r = 5^ + h>^{x) , r,{x) = -la^^ele^^.^^, (2) 




FIG. 1: p' and I are the outgoing fermion and photon and P indicates the incoming fermion. 



where a""^ = 5(7", t'^]- A semicolon and a comma are also used as alternative ways to indicate covariant and partial 
derivatives respectively. We use units h = c~ 1, the signature of rj^^, is —2 , $t = <&s + ^G, 

^s{x)^rj dz^Txiz) ,^G = -\j dz^Yl^xji{z)~-fp^^^[z)]L^^\z) + ]^j dz^j^^k'^, (3) 

L"^{z) = {x" — z°')k^ — {x^ — z^)k" and ^0(2;) satisfies the usual, flat spacetime Dirac equation. 
It is convenient to re-write ([T]) in the form "^{x) — g{x)exp{—ipx)uo{p), where 

9ix) = ^ [7" (Pm + Ki^)Ps. + <^gA^)) + ^] e"^*" , (4) 

and 

^G.f. = -^J dz^{-f^x^i3~-f[ix.^)p'^ + ^lat.p"- (5) 

It also is convenient to focus on the simple process of FigH] We claim that the transition amplitude can be calculated 
by introducing the generalized four-momentum 

P^^ = Pa' + Kp& + ^G,M - ^ (*G - -^^K = Pm + ^M ' (6) 

for the incoming fermion, as (jj]) itself suggests. The part that contains the gravitational field is indicated by P^. 
In ([5]), h", $G,/j and <I>g are quantities that must be calculated, once the metric is known. They are related to the 
Fourier transforms of the corresponding expressions that appear in © , (jH]) and (O . P^ is not on-shell. In fact 



pf^P^ = ml = m'^ + 2 
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where p^^p^ ~ m? because p^ is the momentum of the free fermion represented by ^0(2:) in ©■ The transition 
amplitude is then 

Ml = -iZe7ruo{p')ef,(^)^o9m)MP) , (8) 

where g = p — p' — i, ef^-i represents the polarization of the photon, and Ze is the charge of the fermion. A transition 
amplitude AI2 must be added to Mi to account for the fact that the contraction in (El) is in general accomplished by 
means of g^'^ . It has been repeatedly calculated in the literature and is given by [10| 

A/2 = -iZer-'i\q\)uQip')e^(^^)j,uoip} . (9) 

Ml contains the part p^ of ([6|), that comes from ^Oi and a new part that contains the gravitational contribution due 
to the propagation of the fermion in the field of the source. The total transition amplitude is given by M = Mi -I- M2. 

The calculation now requires that a metric be selected. 

Let us consider the particular instance of a fermion that is propagating with momentum p^ = p, impact parameter 
b > R and X2 = 0, from X3 — — cx) toward a gravitational source of mass M and radius R placed at the origin 
and described by the metric 700 — '^4> tIij = '^4'^ij^ where (/> = — i^. This metric is frequently used in lensing 



problems [ll|,18|. One finds Tq = ~l/2cj)ja°i ,Ti = -1/2(1) j<j'^ and e? = , e? = 1 - (/) , e^ = {l + <j))5{. All 
spin matrices are now expressed in terms of ordinary, constant Dirac matrices. We also assume that the on-shell 
conditions p'p'^ — m? ,£f_i£'^ — remain valid. Extension of the calculation to include different particles, or higher 
order gravitational contributions to p' , £, and ([T]) can be derived to all orders in j^^ . 

The Fourier transforms of the quantities that appear in Q must now be calculated. We obtain 

K{q)Po. = 8n^6iqo)d{q^)6{qy)poGMKoibq,) , /if (gK = S7T^6iqo)S{q^)d{qy)pGMKoibq,) , (10) 

(q) + ^h.M) - , $G,2(<Z) = , $G,3(g) - -87T^6{qo)6{q,)S{qy) (?^ + p) GMKo{bq,) , (11) 
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$G(g) 



-il6TT^d{qo)S{qx)6{qy) 
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Koibqz). 



(12) 



Four- momentum conservation to zeroth order only is required because (flUl) . pT|) and ([T^ are already of 0(7^1,). 
We further approximate the Bessel function Ko{bqz) — \/7r/26(7ze~'"?^[l — 1/86^2 + ...], itself a distribution, by 
Ko{bqz) — S{bqz) and eliminate 6'^{q) from (fTO|) . (fTT|) and ([T2|) . Conservation of energy-momentum will reappear as a 
factor (27r)^(5'*(q) in the expression for the radiated power W defined below. The term $0(9) does however diverge 
even more rapidly than the other Fourier transforms for small momentum transfers and will be dropped in what 
follows. This behavior is well known and is related to the infinite range of the Newtonian potential and the use of a 
plane wave for ^o(a^)- By removing (5^(g) from h" and $G.p(g) we obtain h" and $0./^ of ®- We find 
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We calculate the power radiated as photons in the process of Fig IT] according to the formula [1 



P. 



^-8(^/^^(^-^'-^)^^^^'^^ 



(13) 



(14) 



There are two ways to calculate |Mip. In the first one we replace Pa with Pa in the field- free (7^1^ = 0) expression 
given by E|Mip = Z^e^[— 4m^(p'^p")-|-8(pQp")]. The gravitational contribution to Mi then appears in P^ exclusively. 
We also remove the terms — 32m^(p'^p") -I- 64m^ that do not contain gravitational contributions and therefore refer 
to the kinematically forbidden transition. This yields, to 0(7^^), the expression 



SlMil 



r^2 2 
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^4{p'aPn + »{P^Pn 



(15) 



In a second, alternate approach, we calculate |Afip directly (7^,^ 7^ 0) from ^. By summing over final spins and 
averaging over initial spins and polarizations, we obtain 



S|Mi|2 = -^-^Tr{{p' + m)jp [(p + P + m) ((p + mf + (p + m)P* + H{p + m))] 7''} 



(16) 



where a = 7^*0^ and H — p'^ / p'^ (f){j^ p^ — 7°^"^) -I- {p^/p^j^p'-' — 'y^p^)AGMKi{bqz). On carrying out the traces of 
the Dirac matrices, the contribution from H vanishes. A similar, but simpler calculation, gives |M2p. By further 
eliminating from |Af ^ the terms that refer to the kinematically forbidden transition, we find 



EIMI 



Z^e^ 



-4{p'aP") + 8{paP") 



, 2GM 



-4(p>") -I- Sm^ 



(17) 



The first set of square brackets in P7)) represents the contribution of Mi and coincides with ([T5|) . This supports our 
claim that the generalized momentum P^ introduced in Q leads to the correct value of the transition probability 
by the substitution of p^ with P^ in the field-free expression. The second set of square brackets represents the 
contribution of M2. 

The integration over (Pp' in (|14p is performed by means of the identity J -j^- — J d^p'S{p'^ — m?), while that over 
6 can be carried out by writing the on-shell condition for p' in the form 



M> 



We find 



W ^ 



5(2|P||£| cos 6* -P"P 
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The first term in (|19p represents the contribution of Mi . The radiation spectrum is given by 



dW 2Z2e2 



pz 



(™) 



m^ - £0 {po - p) 



(18) 



(19) 



(20) 



For p > m, p9|) gives 



while it yields 



H^.>™ = ZV^I^, (21) 

b p'^ 



1 ^ lvY^"3 



W,<^ = Z^e^^^('i -'-§], (22) 



for p < m. Equation (J18p and the condition — 1 < cos0 < 1 require that for p > 



m 



fGM\ ( m^\ fGM\ /V A . X 

It follows from (l23l) that the hardest photons are emitted in the forward direction with energy Iq ~ 1Qti{GM /h)p^ / m? 
and power Wpym ^ {\QTiZe)'^{GM/hYp'^ /m? which takes its highest values in the neighborhood of a black hole and 
for high values of p/m. 



For p < m we obtain 



-(T)"'(-^)^'»^-(T)"'(-^)' <-' 



which reduces to the single value £q ^ A'Km{GM /b) when p -^ m. In this case the radiation is still in the forward 
direction and the power radiated 

M/,<„~8(^Ze)2(^)3!!^, (25) 

p'^ 

diverges for small values of p (infrared divergence). This divergence arises as a consequence of the finite energy 
resolution Ae of the outgoing fermion. The process, as calculated, is in fact indistinguishable from that in which 
gravitons with energy < Ae are also emitted and from processes in which vertex corrections are present (virtual 
gravitons emitted and reabsorbed by the external lines of FigdJ. When these additional diagrams are calculated all 
infrared divergences disappear [13[- In the particular case at hand, p in (^5)) is simply replaced by AnGMpo/b, as 
requested by the external field approximation. 

The results (fT9 |) -([22 |) ignore the back-reaction on the background spacetime. We show below that this is not always 
negligible and provide an example of how a very small disturbance in the wave equation can grow rapidly and alter 
the background gravitational field. 

Equation ([T]) requires that '^o{x) be a solution of the field- free Dirac equation and, of course, of the equation 
{ri^^^d^d^ + m^)^o(x) = 0. The approximation procedure still holds true, however, when ^o(a^) satisfies more general 
equations [3|, [ij] . With the addition of a dissipation term, the equation for ^o becomes 

{v'^'d^d^ -f TO^ - 2mado) ^o = , (26) 

where we take a = a|(^o|T'|^o)P = '^{~^^)^ M ^^^ a is a dimensionless, arbitrary parameter, < a < 1, that 
reflects the coupling strength of the dissipation term. When we substitute ^0(2;) = exp{maxo)(t)o{x) into (^5)1 . we 
obtain 

[d^ - d^, + m\l - <j^)] M^) ^ ■ (27) 

An example of a problem with similar behavior is offered by a fluid heated from below. For small temperatures 
gradients the fluid conducts the heat, but as the gradient increases conduction is not sufficient to lead the heat away 
and the fluid starts to convect. In realistic problems the exponential growth of ^0 does not continue indefinitely, but 
is restricted at times xq > t = 1/ma by nonlinearities or dispersive effects that may have been initially neglected. 

The effect of the new solution ^0 on W can be found as follows. We first neglect the change m — > m\/l — a'^ in W 
because in general a < 1. The effect of the exponentially increasing term on (|T4)) then amounts to the transformations 

S{Po-p'q-£o) -^ S{-2ia+Po-p'„-eo) and ^ -^ i3^po-2tay^ - ff3[(pg-4j'2°)2+4pji^2] , where we have used the relation /3 = 

p/pq- M^ has therefore a resonance at po = 2(t of width 4cr^ . Over times xq > t = {am)~^{m/po)~'^{GM/2b)~^GeV~^, 



\If = T^o increases exponentially until the compensating mechanisms mentioned above kick in. For a proton of energy 
Po ^ lOGeV in the field of a canonical neutron star r ~ 3.5 x 10~^^a~^ s. Considerably higher values of r can, of 
course, be obtained for the lighter fermions. As ^ grows, so does the energy momentum tensor associated with it and 
the gravitational field it generates, altering, in the process, the gravitational background. 

Summarizing, external gravitational fields in radiative processes can be included in the calculation of a transition 
probability by simply replacing the momentum p^ of a particle with its generalized version P^ in the corresponding 
expression for the zero-field process. The example given refers to fermions, but the procedure can be extended 
to particles of different spins. An essential point is here that the dispersion relations are altered by the external 
gravitational field and can be calculated if the corresponding wave equations can be solved to 0(7^1^), or higher 
til) ISnSl . The treatment of particle lines in Feynman diagrams therefore necessitates care when external gravitational 
fields are present. It follows, in particular, that kinematically forbidden processes like that of FigH] become physical 
and their transition probabilities can be determined. The calculation of the gravitational contributions are greatly 
simplified and can be extended to higher order in 7^^ . The applications are not confined to fields of a Newtonian type, 
but extend to any gravitational fields. In this respect, the procedure presented goes beyond the results that apply to 
external electromagnetic potentials [9| , not only because the gravitational field has in general ten components rather 
than just four, but also because time-independence is not required in ([H)-®. 

The procedure also yields transition amplitudes and decay rates that, to leading order, are linear in 7^1. and can 
therefore be considerably larger than those normally studied in the literature. These results suggest that some particle 
decay processes in the neighborhood of compact astrophysical objects, or in cosmology, need to be re-evaluated. 

We have also shown that the back-reaction of the fermion on the gravitational background need not be negligible. In 
fact, the addition of a small dissipation term in ((26|) can drastically transform the physical problem over a characteristic 
time r in two respects: i) The power radiated by the fermion acquires a resonant, narrow peak at po '^ 2a and, ii) ^ 
grows exponentially and, in so doing, affects the background gravitational field via its associated energy-momentum 
tensor, until so allowed by the physical circumstances. 
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